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Abstract. Available redundancy among aircraft control surfaces allows for effective wing camber 
modifications. As shown in the past, this fact can be used to improve aircraft performance. To date, 
however, algorithm developments for in-flight camber optimization have been limited. This paper 
presents a perturbational approach for cruise optimization through in-flight camber adaptation. The 
method uses, as a performance index, an indirect measurement of the instantaneous net thrust. As such, 
the actual performance improvement comes from the integrated effects of airframe and engine. The 
algorithm, whose design and robustness properties are discussed, is demonstrated on the NASA 
Dryden B-720 flight simulator. 
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1. INTRODUCTION 

Because of increasing competition among airline 
manufacturers and operators, the world has recently 
seen an all-out effort to reduce airline direct operat- 
ing costs (DOCs). Given that an airline’s net profit is 
a percentage based on the difference between two 
large numbers (revenues and costs), a small drop in 
costs has considerable leverage in an industry with 
annual profits of just about 5 percent (National 
Research Council, 1992). Besides ownership (-50 
percent of DOCs) the second major driver of the 
DOCs of an airline is fuel consumption, which 
accounts for around 18 percent. This shows the 
impact of aircraft performance on an operator’s prof- 
itability, particularly because the operator does not 
usually control production costs. 

In the late 1970’s and 1980’s, many research activi- 
ties were conducted to improve performance aircraft 
as part of NASA programs on aircraft energy effi- 
ciency and advanced fighter technology integration. 
(Gilyard and Espafia (1994) surveys past attempts to 
apply active controls to improve aircraft perfor- 
mance.) The F- 111 mission adaptive wing program 
(Phillips and Smith, 1989) showed the real potential 
of the variable-wing-camber concept to improve air- 
craft performance. Two modes of the F- III mission 
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adaptive wing are of significant interest: (1) the 
cruise camber control mode, designed to use real- 
time adaptive optimization for drag reduction; and 
(2) the maneuver camber control mode, designed to 
maximize lift-drag ratio. The concept of both modes 
was successfully demonstrated in flight, though each 
mode has methodological limitations in practical 
implementation. The maneuver camber control mode 
requires accurate models stored in table lookup form 
onboard the aircraft. The cruise camber control mode 
has an optimization algorithm that is not suited for 
the low levels of drag improvements available on 
transport aircraft. 

More recently, extensive wind-tunnel testing and 
flight experiments with wide-body transport, per- 
formed in Germany at the MBB Company (a mem- 
ber of the European consortium, Airbus Industrie), 
show that continuous camber variations can improve 
the efficiency of even the most advanced wings at 
their best design point (Renken, 1988). The reason is 
that the wing design is the result of a compromise 
between all possible flight conditions. The same 
team shows that coordinated deflections of flaps, 
ailerons, elevators, horizontal stabilizer, and possibly 
leading-edge devices can be used to induce variable 
wing cambering without the penalty of added weight 
or new wing box design (Renken, 1985; Renken, 
1988; Szodruch, 1985). Besides drag improvements 
(potential increase of 3 to 9 percent in lift-drag ratio 



(Renken, 1988)) camber control may also have a 
positive impact on other aspects of the aircraft 
design, such as buffet limit (12-percent increase 
(Renken, 1985; Renken, 1988)) or root bending 
moments (Szodruch, 1985). The multiple benefits of 
the technology have recently prompted the highest 
level of the Airbus consortium to declare that . . 
variable cambering wing features will be mandatory 
for future megatransports” (New Large Aircraft, 
1993). 

The recent introduction of digital fly-by-wire tech- 
nology in modern transport aircraft considerably 
eases the implementation of wing camber optimiza- 
tion techniques. To date, however, algorithm devel- 
opments for in-flight adaptive camber optimization 
have been limited. As such, the development of a 
robust and efficient algorithm to adjust, in flight, the 
aircraft control surfaces configuration is in order. 
This paper considers a perturbational approach for 
cruise optimization through in-flight wing camber 
adaptation. The method is based on the measurement 
of a performance index that includes the engine con- 
tributions. The actual performance improvement 
comes, thus, from the integrated effects of airframe 
and engine. The algorithm is demonstrated on the 
B-720 flight simulator of NASA Dryden Flight 
Research Center. 


symmetric deflection of the outboard ailerons (5 ail ). 
The flight altitude, seen as an optimization con- 
straint, is held by the autopilot, which sends com- 
mands to the elevator (5^,) and the horizontal 
stabilizer (5^). 

From the longitudinal flight equations the horizontal 
acceleration component, a x , is given by 

PI = a = ~[T(FC) cosa- 

* ™ V (1) 

q(v, P)SC D (5„ (/ , FC) - mg siny Jcosy 

where a is the angle of attack, m is the mass of the 
aircraft, q is the dynamic pressure, V 7 is the airspeed, 
p is the air density, S is the wing area, and C D is the 
drag coefficient. FC stands for flying conditions and 
includes all unknown or unmodeled effects of those 
changing conditions, such as weight, center-of- 
gravity position, winds, velocity, altitude, aircraft 
aging, etc. T(FC) corresponds to the engine static 
characteristics, relating thrust (T) with FC. The sur- 
face deflections 5^ and 5 st ^ do not appear explicitly 
in Eq. (1) because these values are implicitly deter- 
mined by the autopilot for a given deflection. 
Under (perfect) cruise conditions (i.e., y = y = 0, 
where y is the path angle) a x = 0, while the total 
derivative, 


Symmetric outboard and inboard ailerons, flaps, ele- 
vator, and horizontal stabilizer are possible surfaces 
available for optimization. The concept is illustrated 
here using symmetric outboard ailerons to produce 
effective wing camber changes together with the 
horizontal stabilizer for pitch-moment compensa- 
tion. This paper addresses the analysis and design 
aspects of the algorithm. Research is currently in 
progress to extend these results to multisurface opti- 
mization. This extension could easily include the 
integration of decision variables from both the air- 
frame and the engine. 
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is equal to zero at the optimal aileron deflection 5*„ 7 
(i.e., for daj&* ai i )ldS ai! = 0), which corresponds to an 
extremum of the acceleration, here used as the PI. 
Notice that the optimization naturally takes into 
account engine thrust changes caused by velocity 
changes. Figure 1 shows the autopilot and the opti- 
mizer loops, where h, h , and h D are, respectively, the 
altitude, its time-derivative, and the engaged altitude. 


2. PROBLEM FORMULATION 

Adaptive optimization techniques with periodic per- 
turbation and direct feedback of a measurable per- 
formance index (PI) allow for direct performance 
optimization without requiring a model (Eveleigh, 
1967; Draper and Li, 1951). In this paper, use is 
made of the acceleration (as a net thrust measure- 
ment) and velocity as indexes. The flight regime 
considered is cruise with fixed power lever angle 
(PLA) position. Similar principles can be applied to 
optimize range or fuel flow (Gilyard and Espana, 
1994). The in-flight optimization of other flight 
regimes, such as uniform climbing, could also be 
envisioned. The algorithm’s viability is demon- 
strated using a single decision variable chosen to be 


PLA (constant) 



Fig. 1. Control diagram with autopilot and optimizer 
loops. 










(P3) The effects of measurement noise on PI are 
negligible unless its power spectrum is concen- 
trated on the excitation frequency co 0 . (This 
may impose conditions on the choice of the 
amplitude A , the frequency G) 0 , or both.) 

Property (P3) is enhanced in practice by using a 
small band pass filter centered on the excitation 
frequency o> 0 . Let cp 0 be the phase angle introduced 
by the cascaded block G f G p (s) at 0) 0 , <p is used to 
compensate for this phase shift. The term 5 0 is the 
trim value at the time of the algorithm engagement. 

For design and analysis purposes, only the informa- 
tion contained in the low-frequency components of 
the signals in the circuit is of interest. For that 
reason, the system in Fig. 2 is transformed to an 
equivalent low-pass network by following the guide- 
lines of reference 8, chapter 9. Figure 3 shows the 
corresponding equivalent network for a quadratic J 
as in Eq. (3). Under each of the blocks of Fig. 3 is 
indicated, in parenthesis, the_ originating block from 
Fig. 2. In the same figure, 5 r represents the (aver- 
aged) low-frequency component of the feedback sig- 
nal 5 r . The parameters of the equivalent system are 
calculated as (Eveleigh, 1967) 


Fig. 3, using standard linear techniques (e.g., root 
locus), one can easily conclude that for K sufficiently 
small and z p > 0, x / > 0 (which is the case in prac- 
tice), the integrator on the feedback path ensures the 
exponential convergence of 5 r to its optimal value 
5* , even under model inaccuracies. The convergence 
to the optimum is, thus, a robust property of the 
algorithm with respect to plant variations. As we 
saw, a small K is also required for the validity of our 
analysis. 

5. APPLICATION TO ACTIVE CAMBER 
OPTIMIZATION 

The algorithm was tested using the NASA Dryden 
B-720 aircraft simulator at a nominal cruise flight 
condition (30,000 ft, 0.8 Mach, and total weight of 
200,000 lb). The frequency to 0 = 0.025 rad/sec is 
purposely made small in comparison with the band- 
width of the aircraft in closed loop with the autopilot. 
With this choice the autopilot ensures a good ride 
quality and diminishes the dynamic effects over the 
optimization constraints. With similar criteria an 
aileron deflection amplitude A = 2° was considered 
acceptable. 
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where |G(yco 0 )| and \j/(0) 0 ) indicate, respectively, 
the modulus and the phase angle of the transfer func- 
tion G at (D 0 . The equivalent network in Fig. 3 is then 
used for the stability analysis and design of the adap- 
tive optimization algorithm. The designer has at his 
disposal the gain K, the phase compensation (p, the 
perturbation signal amplitude A, and the frequency 
co 0 . The transfer functions of the filters G^and G a can 
also be used to refine the design. From the system in 



Fig. 3. Equivalent circuit of block diagram from 
Fig. 2. 


5.1. Desig n 

For the design, the sensor/aircraft dynamics (Fig. 2) 
are neglected by assuming P^ = 1 and = 0. The 
trim point characteristic relating 5 m/ with the PLA 
position (this curve can be determined experimen- 
tally in practice) was fitted with a second-order 
parabola used to calculate the characteristic’s aver- 
age curvature ^PLA' This value r PM is then normal- 
ized by q = l/2pV 2 (281 lb/ft 2 ), giving: V PLA := 
V PlA lq = 7.7E-04. We then calculate the normalized 
curvature for the acceleration used as a PI as 

fac C := ^plaTpla /m =2.85E-05 (11) 

where k PLA = 230 Ib/deg is the thrust-to-throttle gain 
evaluated from the engine data, and m is the total air- 
craft mass. Notice that, given the convergence 
robustness discussed in the previous section, only 
approximated values of the above constants are 
required. 

To keep the design simple, the following first-order 
unitary gain transfer functions were chosen for G f 
and G a : 
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Selecting (O h = to, = (Aq/2 puts to 0 on the band pass of 
G f and out of the band pass of G n . The phase angle 
induced by G f at the chosen ti) 0 is <py = 26.5°. Because 
the dynamics term G p has been neglected, one 




chooses, as shown in Figs. 1 and 2, (p = cp 0 = (P/(the 
phase shift introduced by Gy). From Eq. (10) one has 
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From Eqs. (12) and (13) the open-loop transfer func- 
tion poles of the equivalent system in Fig. 3 are 0, 
(Op/2, and 5cop/2. After normalizing by co 0 , the nor- 
malized open-loop transfer function results in 
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j/co 0 (i/co 0 + 1/2 )(j/co 0 + 5/2) 


? .2 


- G SKqTA 
CO l 4 ®o 


(14) 


5.2. Results 

Figure 3 shows the time histories of a simulated 
speed maximization flight experiment at 30,000 ft 
and 0.8 Mach (approximately 796 ft/sec). The auto- 
pilot regulates the altitude while the PLA is kept at 
62° corresponding, approximately, to an engine 
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where G is the combined gain of the blocks ( J ), (Gy), 
(M), G fl , and KJs of Fig. 3 put in cascade. The nor- 
malized dynamics of the performance optimization 
algorithm is characterized by the closed-loop poles 
of the transfer function (14). Those poles can be con- 
veniently placed by an appropriate choice of the nor- 
malized gain G. The critical normalized gain for a 
deadbeat response can be shown to be G r = 0.1408, 
with a corresponding coincident dominant (and nor- 
malized) closed-loop poles at -0.24 and a single 
normalized fast pole at -2.52. For the case under 
consideration with the average curvature T given by 
Eq. (1 1) and the selected A and co Q , the design results 
in a pair of dominant poles at -0.006 rad/sec and a 
fast pole at -0.063 rad/sec. The algorithm gain, K = 
0.196 is calculated using Eq. (14) as a function of 
the selected G and the rest of the design parameters. 
This completes the algorithm design. Equation (14) 
also gives an easy way to update K with the dynamic 
pressure q . Notice that for a given design criterion 
(in this paper we chose a critical design), the adapta- 
tion gain, the corresponding closed-loop poles, and, 
consequently, the algorithm convergence speed are 
all proportional to the excitation frequency co 0 . 

Given the small accelerations involved, the airspeed 
velocity may be a more suitable performance index 
than the acceleration itself. For this case the corre- 
sponding adaptation gain K vd is calculated using the 
fact that the co 0 -component of the measured signal V x 
satisfies 


dV* 

dt 


Vcosy = 


a x 
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Thus, to use V x as the measured PI, K has to be 
modified by a factor of JcoJ and the phase (p delayed 
90°. This gives K vtl = O.196co 0 and cp = 26.5° - 90° = 
-63.5°. 


Velocity, ft/sec 



Altitude, ft 
30.02 x 10 3 



Normal acceleration, g 



Fig. 4. Time histories from B-720 velocity optimiza- 
tion; Mach 0.8 at 30,000 ft; 200,000-lb gross weight. 



thrust level of 15,000 lb. The algorithm is engaged at 
r = 120 sec. Figure 4 displays the trim command 
(8 nit ) sent to the aileron, the excitation signal super- 
imposed on it, and the corresponding horizontal sta- 
bilizer deflection (S irff ) produced by the autopilot. 
Figure 4 also shows the actual speed time history 
indicating a speed increase of approximately 20 
ft/sec. As seen in the figure, the periodic perturbation 
induces an altitude oscillation with the same funda- 
mental period and an amplitude ranging from 10 to 
20 ft. Finally, the control law is shown to have a min- 
imal effect on the normal acceleration (less than 
5E-03 g in absolute value) with negligible conse- 
quences on the ride quality. 

6. CONCLUDING REMARKS 

The possibility of continuously improving aircraft 
performance using conventional control surfaces for 
wing camber adaptation is demonstrated for a simu- 
lated transport aircraft. To take advantage of the full 
potential of the variable-camber concept, an algo- 
rithm based on a direct measurement of performance 
is required. In this paper a periodic perturbation 
technique is explored and its viability demonstrated. 
Besides its ease of implementation, the main advan- 
tages of the proposed technique are (1) its noise 
insensitivity and (2) its robustness with respect to 
measurement biases and model inaccuracies. Fur- 
thermore, the algorithm’s design is simple and 
straightforward. Applying the algorithm to the 
velocity optimization problem results in a significant 
increase in cruise velocity for a simulated B-720 
transport aircraft with in-flight adaptive reposition- 
ing of aileron/horizontal stabilizer deflections. The 
same algorithm can also optimize other performance 
indices such as range, acceleration (net thrust) or 
fuel flow. By using the proposed algorithm, perfor- 
mance improvement can be achieved for different 
flight conditions, weight changes, aircraft aging, and 
center-of-gravity position changes. The principle can 
be extended to multivariable optimization, which 
makes the algorithm suitable for application to air- 
frame-propulsion integrated optimization. The fact 


that future aircraft will incorporate active control 
features in their aerodynamics and propulsion con- 
trol effectors further sustains the significance of the 
results presented in this paper. 
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